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Abstract. We report on results from about 30 hours of livetime with the Goldstone 
Lunar Ultra- high energy neutrino Experiment (GLUE). The experiment searches for 
I < 10 ns microwave pulses from the lunar regolith, appearing in coincidence at two large 

■T^lj" . radio telescopes separated by about 22 km and linked by optical fiber. The pulses can 

CN \ arise from subsurface electromagnetic cascades induced by interactions of up-coming 

. ~ 100 EeV neutrinos in the lunar regolith. A new triggering method implemented 

after the first 12 hours of livetime has significantly reduced the terrestrial interference 
background, and we now operate at the thermal noise level. No strong candidates are 
yet seen. We report on limits implied by this non-detection, based on new Monte Carlo 
estimates of the efficiency. Wc also report on preliminary analysis of smaller pulses, 
O \ where some indications of non-statistical excess may be present. 

to 



I INTRODUCTION 

Recent accelerator results [1,2] have confirmed the 1962 prediction of Askaryan 
[3,4] that electromagnetic cascades in dense media should produce strong coher- 
ent pulses of microwave Cherenkov radiation. These confirmations strengthen the 
motivation to use this effect to search for cascades induced by predicted diffuse 
backgrounds of high energy neutrinos, which are associated with the presence of 
> 10^° eV cosmic rays in many models. At neutrino energies of about 100 EeV (1 
EeV = 10^^ eV), cascades in the upper 10 m of the radio-transparent lunar regolith 
result in pulses that are detectable by large radio telescopes at earth [5,6]. One 
prior experiment has been reported, using the Parkes 64 m telescope [7] with about 
10 hours of livetime. 

At frequencies above 2 GHz, ionospheric delay smearing is unimportant, and 
the signal should appear as highly linearly-polarized, band-limited electromagnetic 
impulses [8-10]. However, since there are many anthropogenic sources of impul- 



sive radio emission, the primary problem in detecting such pulses is eliminating 
sensitivity to such interference. 

Since 1999 we have been conducting a series of experiments to establish tech- 
niques to measure such pulses, using the JPL/NASA Deep Space Network antennas 
at Goldstone Tracking Facility near Barstow, California [11]. We employ the 70 m 
and 34 m telecommunication antennas (designated DSS14 and DSS13 respectively) 
in a coincidence-type system to solve the problem of terrestrial interference, and 
this approach has proven very effective. Since mid-2000, the project has moved 
into a new status as an ongoing experiment, and receives more regularly scheduled 
observations, subject to the constraints imposed by the spacecraft telecommunica- 
tions priorities of the Goldstone facility. 

Although the total livetime accumulated in such an experiment is a relatively 
small fraction of what is possible with a dedicated system, the volume of material 
to which we are sensitive, a significant fraction of the Moon's surface to ~ 10 m 
depth, is enormous, exceeding 100,000 km^ at the highest energies. The result- 
ing sensitivity is enough to begin to constrain some models for diffuse neutrino 
backgrounds at energies near and beyond 10^*^ eV. We report on the status of the 
experiment, and astrophysical constraints imposed by limits from about 30 hours 
of livetime. We are also improving our understanding of the emission geometry and 
detection sensitivity through simulations, and describe initial results in extending 
our sensitivity to pulses of lower amplitude. 

II DESCRIPTION OF EXPERIMENT 

The lunar regolith is an aggregate layer of fine particles and small rocks, thought 
to be the accumulated ejecta of meteor impacts with the lunar surface. It consists 
mostly of silicates and related minerals, with meteoritic iron and titanium com- 
pounds at an average level of several per cent, and traces of meteoritic carbon. 
It has a typical depth range of 10 to 20 m in the maria and valleys, but may be 
hundreds of meters deep in portions of the highlands [12]. It has a mean dielectric 
constant of e ~ 3 and a density of p ^ 1.7 gm cm~^, both increasing slowly with 
depth. Measured values for the loss tangent vary widely depending on iron and 
titanium content, but a mean value at high frequencies is tan 5 ~ 0.003, implying 
a field attenuation length at 2 GHz of {a)~^ = 9 m [13]. 

A Emission geometry Signal Characteristics 

In Fig. 1 we illustrate the signal emission geometry. At 100 EeV the interaction 
length Lint of an electron or muon neutrino for the dominant deep inelastic hadronic 
scattering interactions (averaging over the charged and neutral current processes) 
is about 60 km [15] (i?^ = 1740 km). Upon interaction, a ~ 10 m long cascade 
then forms as the secondary particles multiply, and compton scattering, positron 
annihilation, and other scattering processes then lead to a ~ 20% negative charge 



excess which radiates a cone of coherent Cherenkov emission at an angle of 56°, with 
a FWHM of 1°. The radiation propagates in the form of a sub-ns pulse through 
the regolith to the surface where it is refracted upon transmission. 

Because the angle for total internal reflection (TIR) of the radiation emitted from 
the cascade is to first order the complement of the Cherenkov angle, we consider for 
the moment only neutrinos that cascade upon emerging from a penetrating chord 
through the lunar limb. Under these conditions the typical neutrino cascade has 
an upcoming angle with respect to the local surface of 



which implies a mean of 6up ~ 1° at 10^'^ eV. 

At the regolith surface the resulting microwave Cherenkov radiation is refracted 
strongly into the forward direction. Scattering from surface irregularities and de- 
magnification from the interface refraction gradient fills in the Cherenkov cone, 
and results in a larger effective area of the lunar surface over which events can be 
detected, as well as a greater acceptance solid angle. These effects are discussed in 
more detail in section III. A. 



The antennas employed in our search are the shaped-Cassegrainian 70 m antenna 
DSS14, and the beam waveguide 34 m antenna DSS13, both part of the NASA 
Goldstone Deep Space Network (DSN) Tracking Station. DSS13 is located about 
22 km to the SSE of DSS14. The S-band (2.2 GHZ) right-circular-polarization 
(RCP) signal from DSS13 is filtered to 150 MHz BW, then downconverted with an 
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FIGURE 1. Schematic of the geometry for lunar neutrino cascade event detection. 



intermediate frequency (IF) near 300 MHz. The band is then further subdivided 
into high and low frequency halves of 75 MHz each, and no overlap. These IF 

signals are then sent via an analog fiber-optic link to DSS14. At DSS14, the dual 
polarization S-band signals are downconverted with the same 300 MHz IF, and 
bandwidths of 100 MHz (RCP) and 40 MHZ (LCP) are used. A third signal is also 
employed at DSS14: a 1.8 GHz (L-band) feed which is off-pointed by ~ 0.5° is used 
as a monitor of terrestrial interference signals; the signal is downconverted in the 
same manner as the other signals and has a 40 MHz bandwidth. 

C Trigger system 

The experimental approach in our initial 12 hours of observations was to use a 
single antenna trigger with dual antenna data recording [1 1] . This was accomplished 
by using the local S-band signals as DSS14 to form a 2-fold coincidence with an 
active veto from the L-band interference monitor. Since any system with an active 
veto is subject to potential unforeseen impact on the trigger efficiency, we have 
now developed an approach which utilizes signals from both antennas to form a 
real-time dual-antenna trigger, with no active veto. 

Fig. 2 shows the layout of the trigger. The four triggering signals from the two 
antennas are converted to unipolar pulses using tunnel-diode square-law detectors. 
Stanford Research Systems SR400 discriminators are used for the initial threshold 
level, and these are set to maintain a roughly constant singles rate, typically 0.5-1 
kHz/chan for DSS14 and 30 kHz/chan for DSS13 (DSS13's rate is higher due to a 
lower threshold, compensating for the reduced aperture size). A local coincidence 
is then formed for each antenna's signals. The DSS14 coincidence between both 
circular polarizations ensures that the signals are highly linearly polarized, and the 
DSS13 coincidence helps to ensure that the signal is broadband. 

Fig. 3 indicates the timing sequence for a trigger to form (negative logic levels 
are used here). A local coincidence at DSS14, typically with a 25 ns gate, initiates 
the trigger sequence. After a 65 /is delay, a 150 /is gate is opened (the delays 
compensate for the 136 fxs fiber delay between the two antennas). This large time 
window encompasses the possible geometric delay range for the moon throughout 
the year. Use of a smaller window is possible but would require delay tracking and 
a thus more stringent need for testing and reliability; use of a large window avoids 
this and a tighter coincidence can then be required offline. 

If a 25 ns local coincidence now forms between the two DSS13 signals within the 
allowed 150 /xs window, a trigger is formed. The sampling scopes are then triggered, 
and a 250 record, sampled at 1 Gs/s, is stored. The average trigger rate, due 
primarily to random coincidences of thermal noise fluctuations, is about 1.6 mHz, 
or 1 trigger every 5 minutes or so. Terrestrial interference triggers are uncommon 
(a few percent of the total), but can occasionally increase in number when a large 
burst of interference occurs at either antenna, with DSS14 more sensitive to this 
effect. The deadtime per event is about 6 s; thus on average we maintain about 



99% livetime during a run. 



Ill ESTIMATED SENSITIVITY 

Estimates of the sensitivity of radio telescope observations usually involve sys- 
tems that integrate total power for some time constant At which is in general much 
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FIGURE 2. The GLUE trigger system used for the lunar neutrino search. 
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FIGURE 3. A timing diagram for the GLUE trigger system. 



longer than the antenna's single temporal mode duration which is given by the in- 
verse of the bandwidth: r = (Ai/)~^. Since the pulses of interest in our experiment 

arc much shorter than this time scale, the observed pulse structure of induced 
voltage in the antenna receiver is determined only by the bandpass function; that 
is, the pulses are band-limited. Thus the typical dependence of sensitivity on the 
factor V AtAu does not obtain; this factor is always unity in band-limited pulse 
detection. 

Because much of the theoretical work in describing such pulses has been done in 
terms of field strength rather than power, we analyze our sensitivity in these terms 
as well. Such analysis is also compatible with the receiving system, which records 
antenna voltages proportional to the incident electric field, and leads to a more 
linear analysis. It also yields signal-to-noise ratio estimates which are consistent 
with Gaussian statistics, since thermal noise voltages are described by a Gaussian 
random process. 

The expected field strength per unit bandwidth from a cascade of total energy 
Wt can be expressed as [8-10]: 

Eo V m"^ MHz-^) = — — , , , 2 

where R is the distance to the source in m, u is the radio frequency, and the 
decohcrencc frequency is z/q ~ 2500 MHz for regolith material (z/q scales mainly 
by radiation length). For typical parameters in our experiment, a 10^^ eV cascade 
will result in a peak field strength at earth of ~ 0.5 //V m"^ for a 70 MHz BW. 
Equation 2 has now been verified to within factors of 2 through accelerator tests [1,2] 
using silica sand targets and 7-ray-bunch-induced cascades with Wt < 10^^ eV per 
bunch. 

Given that the use of a dual antenna trigger has virtually eliminated the problem 
of terrestrial interference that was the primary limitation to the sensitivity of the 

one previous experiment [7], we can now express the minimum detectable field 
strength Emin for each antenna in terms of the induced signal and the thermal 
noise background. 

The expected signal strength Eq induces a voltage at the antenna receiver given 

by 

Vs = heEoAi^ (3) 
where the antenna effective height he is given by [14] 



he = 2^1^ cose, (4) 



where Za is the antenna radiation resistance, 77 and A are the antenna efficiency 
and area, respectively, Zq = 377 VL is the impedance of free space, and 9, the 



polarization angle of the antenna with respect to the plane of polarization of the 
radiation. 

The average thermal noise voltage in the system is given by 



AkT^ysZT^u . (5) 

Here k is Boltzmann's constant, Tgys is the system thermal noise temperature, and 
Zt the termination impedance of the receiver. If we assume that Z^ ~ Zt then the 
resulting SNR S is 



kTgygZ( 

The minimum detectable field strength (Eq — > Emin) is then given by 



S = ^ ^ E,cose,J-^—— . (6) 



V rjAAu cos Op 

Combining this with equation 2 above, the threshold energy for pulse detection is 



WthriEeV) ~ 4.0 
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(8) 



For the lunar observations on the limb, which make up about 85% of the data 
reported here, Tgys — 110 K, i/ = 2.2 GHz, and the average Ai> ~ 70 MHz. For 
the 70 m antenna, with efficiency rj ~ 0.8, the minimum detectable field strength 
is E^in ~ 1.2 X 10-^ V m"^ MHz"^ for cos6'p = 0.7. The estimated threshold 
energy for these parameters is Wthr = 2.8 x 10^^ eV, assuming a detection level 
of = 5 (5(t) per IF at DSS14 (with a somewhat lower requirement at DSS13 in 
coincidence) . 



A Monte Carlo results 

To estimate the effective volume and acceptance solid angle as a function of 
incoming neutrino energy, events were generated at discrete neutrino energies, in- 
cluding the current best estimates of both charged and neutral current cross sec- 
tions [15], and the Bjorken-y distribution. Both electron and muon neutrino in- 
teractions were included, and Landau-Pomeranchuk-Migdal effects in the shower 
formation were estimated [10]. At each neutrino energy, a distribution of cascade 
angles and depths with respect to the local surface was obtained, and a refraction 
propagation of the predicted Cherenkov angular distribution was made through the 
regolith surface, including absorption and reflection losses and a flrst order rough- 
ness model. Antenna thermal noise fluctuations were included in the detection 
process. 
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FIGURE 4. The microwave Cherenkov radiation pattern from an event in the lunar regoHth. 

A portion of the simulation is shown in Fig. 4. Here the flux density is shown as 

it would appear projected on the sky, with (0,0) corresponding to the tangent to 
the lunar surface in the direction of the original cascade. The units are Jy (1 Jy = 
10^^^ W m^^ Hz^^) as measured at earth, and the plot is an average over several 
hundred events at different depths and a range of 9 up consistent with a 10^° eV 
neutrino interaction, averaging over inelasticity effects and a mixture of electron 
and muon neutrinos consistent with decays from a hadronic tt^ source. 

Although the averaging has broadened the distribution somewhat, a typical cas- 
cade still produces a flux density pattern of comparable angular size. The angular 
width of the pattern directly increases the acceptance sohd angle, and the angular 
height increases the annular band of the lunar surface over which neutrino events 
can be detected, as indicated in Fig. 1. The net effect is that, although the specific 
flux density of the events are lowered somewhat by refraction and scattering, the 
effective volume and acceptance solid angle are significantly increased. The neu- 
trino acceptance solid angle, in particular, is about a factor of 50 larger than the 
apparent solid angle of the moon itself. 

B EHE cosmic rays 

We have noted above that the refraction geometry of the regolith favors emission 
from cascades that are upcoming relative to the local regohth surface. Thus to first 
order EHE cosmic ray events, which cascade within a few tens of cm as they enter 
the regolith, will not produce detectable pulses since their emission with be totally 
internally refiected within the regolith. This effect has been now demonstrated in 
an accelerator experiment [2]. 

This conclusion does not account for several effects however. These effects are 
illustrated in Fig. 5. In Fig. 5A, the varying surface angles due surface roughness on 
scales greater than a wavelength will lead to escape of some radiation from cosmic 



ray cascades. Fig. 5B illustrates that, because the Cherenkov angular distribution 
is not infinitely narrow arround the Cherenkov angle (FWHM ~ 1°), emission at 
angles larger then the Cherenkov angle can escape total internal reflectance. In 
Fig. 5C cascades from cosmic rays that enter along ridgelines can encounter a 
change in slope of the local surface that results in more efficient transmission of the 
radiation. 

Even if total internal reflection strongly suppresses detection of cosmic rays in 
cases A and B, the latter case C of favorable surface geometry along ridgelines or 
hilltops will lead to some background of EHE cosmic ray events. We have not as 
yet made estimates of this background.^ 

Fig. 5D shows the formation zone aspect of the process of Cherenkov emission 
from near-surface cosmic ray cascades. This constraint may suppress Cherenkov 
production even if surface roughness and the width of the Cherenkov distribution 
would otherwise favor some escape of emission. 
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FIGURE 5. Various effects associated with EHE cosmic ray hadron interactions. 

It has now been conclusively shown [16] that coherent Cherenkov emission is a 
process involving the bulk dielectric properties of the radiating material. Cherenkov 
radiation is induced over a macroscopic region of the dielectric (with respect to the 
scale of a wavelength), and does not even require that the charged particles enter 
the dielectric for radiation to be produced — a proximity of several wavelengths or 
less is sufficient [17]. A corollary to this result is that a cascade travelling along 
very near a boundary of the dielectric will not radiate (or radiate only weakly) into 
the hemisphere with the boundary. 

These conclusions also apply of course to the fraction of neutrinos which interact on entering 
the regolith as well as those which interact near their projected exit point. Thus we have not yet 
accounted for all of the possible neutrino events as well as the cosmic ray background events. 



Thus in the case of a cosmic ray entering the regohth at near grazing incidence 
(say within ~ 1°) the resulting cascade reaches maximum within ~ 3 cm of the 
surface, still less than a wavelength for S-band observations. We therefore expect 
that the suppression of Cherenkov emission in such events significantly reduces 
our sensitivity to cosmic rays. Such effects have not been included yet in other 
estimates [9,18] of the cosmic ray detection efficieny of such experiments. 

IV RESULTS 
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FIGURE 6. Plot of model neutrino fluxes and limits from the Fly's Eye experiment and the 
present work. 

Figure 6 plots the predicted fluxes of EHE neutrinos from a number of models 
including AGN production [23] gamma-ray bursts [19], EHE cosmic-ray interac- 
tions [22], topological defects [25,21], and the Zq burst scenario [24]. Also plotted 
are limits from about 70 days of Fly's Eye livetime [20] (accumulated in several 
years of runtime), which apply only to electron neutrino events. 

Our initial 90% CL limit, for 30 hours of livetime is shown plotted with diamonds 
(see also Table 1), based on the observation of no events above an equivalent 5a 
level amplitude (referenced to the 70 m antenna) consistent with the direction of 
the moon. These limits assume a monoenergetic signal at each energy; thus they 
are differential limits and independent of source spectral model, and represent the 
most conservative limits we can apply. Our limits just begin to constrain the highest 
topological defect model [25] for which we expected a total of order 1-2 events. 



TABLE 1. Differential limits on mono-energetic EHE neutrino fluxes. 



Energy (eV) lO^'^ 3 x lO^^ 10^° 3 x 10^" 10^^ 3 x 10^^ 10^^ 3 x 10^^ lO^^ 

\ogio{E^dF/dE) -3.14 -3.66 -3.92 -3.91 -3.73 -3.42 -3.03 -2.66 -2.30 
(GeV cm~^ s~^ sr~^) 
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FIGURE 7. Low amplitude event histogram of GLUE data in delay with respect to the moon 
center delay. The overlain histogram is the expected background level. An excess occurs near 
zero delay, with a 2 microsecond offset. 

In addition to the limits set above from the non-observation of events above, we 
have also analyzed events which triggered the system, but did not pass our more 
stringent software amplitude cuts. A sample of events was prepared by applying 
our standard cuts to remove terrestrial interference events. We required somewhat 
tighter timing that the hardware trigger, as well as band-limited pulse shape, but 
allowed smaller amplitudes, typically corresponding to ~ A.ba at DSS14, and about 
3(T at DSS13. The results are shown in Fig. 7, where the passing events have been 
binned according to their delay timing with respect to the expected delay from 
an event at the center of the moon. The background level (solid line) has been 
determined by randomizing the UT of the events and indicates the somewhat non- 
uniform seasonal coverage of our observations. 

An excess is observed in the vicinity of zero delay where the lunar events are 
expected to cluster. At present there is a ~ 2/is offset from zero delay; this is 
too large to be accounted for by differential delays to the lunar limb, which can 
produce offsets of several hundred ns. Further study of the low amplitude events 
is in progress. 




V CONCLUSIONS 



We have developed a robust system for observing microwave pulses produced 
in the lunar regolith by electromagnetic particle cascades above ~ 10^^ eV. We 
have operated this system to achieve a livetime of 30 hours, with no large apparent 
signals detected to date. We have set conservative upper limits on the diffuse 
cosmic neutrino fluxes over the energy range from 10^^~^^ eV. We have also begun 
to analyze smaller events and have some preliminary indications that a signal may 
be present, but requiring further study. 
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